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We investigate the intergenerational shape dynamics of single Caulobacter crescentus cells using 
a novel combination of imaging techniques and theoretical modeling. We determine the dynam¬ 
ics of cell pole-to-pole lengths, cross-sectional widths, and medial curvatures from high accuracy 
measurements of cell contours. Moreover, these shape parameters are determined for over 250 cells 
across approximately 10000 total generations, which affords high statistical precision. Our data and 
model show that constriction is initiated early in the cell cycle and that its dynamics are controlled 
by the time scale of exponential longitudinal growth. Based on our extensive and detailed growth 
and contour data, we develop a minimal mechanical model that quantitatively accounts for the cell 
shape dynamics and suggests that the asymmetric location of the division plane reflects the distinct 
mechanical properties of the stalked and swarmer poles. Furthermore, we find that the asymmetry 
in the division plane location is inherited from the previous generation. We interpret these results 
in terms of the current molecular understanding of shape, growth, and division of C. crescentus. 


Cell shape both reflects jT) and regulates [2] biological 
function. The importance of cell shape is exemplified by 
bacteria, which rely on specific localization of structural 
proteins for spatiotemporal organization pf). Bacteria 
take forms resembling spheres, spirals, rods, and cres¬ 
cents. These shapes are defined by cell walls [1] consist¬ 
ing of networks of glycan strands cross-linked by peptide 
chains to form a thin peptidoglycan meshwork [5] . Super¬ 
resolution imaging is now revealing the internal positions 
of associated proteins [9] . These include cytoskeletal pro¬ 
teins such as MreB, a homolog of actin EM! , interme¬ 
diate filament-like bundles of CreS (crescentin) [TL, T2] . 
and FtsZ, a homolog of tubulin m- However, due to 
the inherently stochastic nature of molecular processes, 
understanding how these proteins act collectively to ex¬ 
ert mechanical stresses and modulate the effects of turgor 
pressure and other environmental factors requires com¬ 
plementary methods such as high-throughput, quantita¬ 
tive optical imaging. 

Multigenerational imaging data for bacterial cells can 
now be obtained from microfluidic devices of various de¬ 
signs [T4j-|T8] . Still, a common limitation of most devices 
is that the environmental conditions change throughout 
the course of the experiment, particularly as geometric 
growth of the population results in crowding of the ex¬ 
perimental imaging spaces. We previously addressed this 
issue by engineering a C. crescentus strain in which cell 
adhesion is switched on and off by a small molecule (and 
inducible promoter) Il9j . allowing measurements to be 
made in a simple microfluidic device [TM22] . This tech¬ 
nology allows imaging >100 generations of growth of an 
identical set of 250-500 single cells distributed over ~25 
fields of view. Thus cell density is low and remains con¬ 
stant. These studies afforded sufficient statistical preci¬ 
sion to show that single C. crescentus cells grow exponen¬ 


tially in size and divide upon reaching a critical multiple 
(«1.8) of their initial sizes jT5]. Satisfaction of a series of 
scaling laws predicted by a simple stochastic model for 
exponential growth indicates that these dynamics can be 
characterized by a single time scale [19J 123] . 

In this paper, we use more advanced image analysis 
methods to extract cell shape contours from these data. 
The resulting geometric parameters, together with math¬ 
ematical models, provide insights into growth and divi¬ 
sion in C. crescentus and the plausible role of cell wall 
mechanics and dynamics in these processes. Specifically, 
we identify natural variables for tracking cell dynamics, 
and develop a minimal mechanical model that shows how 
longitudinal growth can arise from an isotropic pressure. 
We then examine the dynamics of cell constriction and 
unexpectedly find that it is governed by the same time 
constant as exponential growth. This important find¬ 
ing can be understood in terms of an intuitive geometric 
model that relates the constriction dynamics to the kinet¬ 
ics of the growth of septal cell wall. We further suggest 
that the site of constriction can arise from differences in 
materials properties of the poles and show that it is es¬ 
tablished in the previous generation- i.e., the location of 
the site of division can be predicted before formation of 
the divisome. We relate our results to the known dynam¬ 
ics of contributing molecular factors and existing models 
for bacterial growth and division. 


RESULTS 

The length is sufficient to characterize the expo¬ 
nential growth of each cell. Various techniques have 
been put forth to analyze cell morphology gathered from 
single cell images [23]. Recent work on image analysis 
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FIG. 1. Determination of cell contour and definition of shape parameters, (a) A representative phase contrast image 
of one field of view. The solution flow in the microfluidic channel is from bottom to top. (b) Zoomed image of the yellow 
highlighted cell from a and its splined contour, (c) Schematic of a contour illustrating the shape parameters. The cell medial 
axis is calculated from pole to pole; it defines both cell length l(cf)) and radius of curvature which lead directly to the 

spanning angle #(</>)• The cell width u) is a parametric quantity, calculated as the length of the rib perpendicular to the 
medial axis at a specified distance from the stalked pole, The location of the global minimum of the width Wmin (<A) 

(purple line) can be used to segment the cell into stalked ( st , red) and swarmer ( sw , blue) portions. 


of single cells has attempted to optimize two problems: 
separation of distinct (but potentially overlapping) cells 
and accurate determination of the edge of each cell [25] , 
Because crowding is not an issue in our setup, we could 
focus solely on constructing an algorithm to delineate 
each cell contour accurately and precisely. As shown in 
Fig. la,b and described in the Methods section, we first 
segment each cell using pixel-based edge detection simi¬ 
lar to |26l , then perform spline interpolation to determine 
the cell contour at sub-pixel resolution. The sequence of 
such images for each single cell constitutes a trajectory in 
time t that serves as the basis for quantitative analysis. 
Division events are then detected in an automated fash¬ 
ion using custom Python code, and used to divide time 
trajectories for each cell into individual generations. 

All data shown here were obtained by observing 260 
single C. crescentus cells perfused in complex medium 
(peptone-yeast extract; PYE) at 31°C over the course 
of 2 days (corresponding to 9672 separate generations). 
Under these conditions, the mean population growth rate 
and division time remain constant, so we treat the tra¬ 
jectories of individual generations as members of a single 
ensemble. In other words, we segment each cell trajec¬ 
tory by generation and take the resulting initial frame 
(i.e., immediately following division) as t = 0 minutes. 
In order to average over the ensemble, we then bin quan¬ 
titative information according to time since division, t, 
normalized by the respective division time r. The nor¬ 
malized time, <f> = t/r, serves as a cell-cycle phase vari¬ 
able. 

For our quantitative analysis, we focus on a set of three 
intuitive and independent parameters that characterize 
cell shape at each stage of growth: length l, width w, 


and radius of curvature R (Fig. lc). They are calculated 
directly from each splined contour as follows (see also 
Supplementary Fig. 1): 

• We define the length, /(</>), as the pole-to-pole dis¬ 
tance along the contour of the cell medial axis at 
the normalized time </> (Fig. 2a). 

• We assign a single radius of curvature, R(4>), to 
each cell based upon the best-fit circle to the me¬ 
dial axis (Fig. 2b). Although stalked (i? st (</>)) and 
swarmer (R sw (</>)) portions may be described by 
different radii of curvature toward the end of the 
cell cycle, the average radius obtained by averaging 
the contributions of each portion yields the same 
value, i.e., ((I? st ((/>) + R sw ((p))/ 2) ~ (i?(<£)) (see 
Supplementary Fig. 2c). 

• We define the width, w(cf>,u ) as the length of the 
perpendicular segment spanning from one side of 
the cell contour to the other at each position u(<p) 
along the medial axis, which runs from u = 0 at 
the stalked pole to it = Z at the swarmer pole. 
Furthermore, we spatially averaged the width over 
positions along the medial axis, w{</>), to obtain a 
characteristic width at each time point (Fig. 2c). 

The mean division time is (r) = 73± 7 min, where (...) 
indicates a population average. We find that (Z(</>)) in¬ 
creases exponentially with time constant ( k ) _1 = 125 ±8 
min, essentially the same time constant that we previ¬ 
ously observed for the cross-sectional area m, while 
(w(</))) and (!?(</>)) remain approximately constant for 
0 < cj) < 0.5 and each shows a dip for 0.5 < (j) < 0.9 
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when cell constriction becomes prominent. The sharp 
rise in (R(if))) seen for <j> > 0.9 results from independent 
alignment of the stalked and swarmer portions with the 
microfluidic flow as they become able to move indepen¬ 
dently (i.e., fluctuate easily about the plane of constric¬ 
tion). These observations confirm the assumptions in 
m that the length is sufficient to describe the growth 
of cell size. Moreover, we can track the dynamics of the 
spanning angle, 6 : using the relation l = R6. 

Mechanical model for cell shape and growth. 
There are many details of cell growth and shape that 
require interpretation. For example, it is not obvious 
a priori that growth should be almost exclusively lon¬ 
gitudinal. Therefore, we have developed a minimal me¬ 
chanical model that can explain these observations. We 
parametrize the geometry of the cell wall by a collection 
of shape variables {qi(t)}, where q\ = (R), q 2 = (w), and 
53 = (6) are the parameters introduced above (Fig. lc). 
As the cell grows in overall size, we postulate that the rate 
of growth in the shape parameter qi(t) is proportional to 
the net decrease in cell wall energy, -E({g»(t)}), per unit 
change in q^t) [2ZI ;M|- Assuming linear response, the 
configurational rate of strain, ^, is proportional to 

the corresponding driving force F t = — dE/dqt , in anal¬ 
ogy with the constitutive law of Newtonian flow [25] : 


1 dqj 

qi dt 




(1) 


where the constant <f>j describes the rate of irreversible 
flow corresponding to the variable qt(t). According to 
equation Q, exponential growth occurs if F, is constant, 
whereas (p(t) reaches a steady-state value if Fi(qi) = 0 
along with the condition dFi/dqi < 0. It thus remains to 
specify the form of E. 

For a C. crescentus cell of total volume V and surface 
area A, our model for the total energy in the cell wall is 
given by 


E(R,W,9) — —PV + J dA 7 + i? w idth+A'cres+-E'div ) (2) 

where P is a constant pressure driving cell wall expan¬ 
sion; 7 is the tension on the surface of the cell wall; -E w idth 
is the energy required to maintain the cell width; E cres 
represents the mechanical energy required to maintain 
the crescent cell shape; is the energy driving cell wall 
constriction. Traditionally P was taken to be the turgor 
pressure E3; while the importance of the turgor pressure 
has recently been questioned [3D], an effective pressure 
must still arise from the synthesis and insertion of pep- 
tidoglycan strands that constitute the cell wall. We note 
that a purely elastic description of cell wall mechanics 
would lead to a curvature-dependent surface tension m- 
However, if growth is similar to plastic deformation, the 
tension is uniform [3]. The effective tension in our model 


depends on the local surface curvatures through the en¬ 
ergy terms -E w idth and E cies , that describe harmonic wells 
around preferred values of surface curvatures. 

The mechanical energy for maintaining width is given 

by 


Awi d t h 2 J dA {w/2 R m ) ’ (3) 

where the constant R m is the preferred radius of cur¬ 
vature, k m is the bending rigidity and dA is a differen¬ 
tial area element 33]. Contributions to k m can come 
from the peptidoglycan cell wall as well as membrane- 
associated cytoskeletal proteins like MreB, MreC, RodZ, 
etc., which are known to control cell width [ZHS] - 

In addition to maintaining a constant average width, 
C. crescentus cells exhibit a characteristic crescent shape, 
which relies on expression of the intermediate filament- 
like protein crescentin m ■ Although the mechanism by 
which crescentin acts is not known, various models have 
been proposed, including modulation of elongation rates 
across the cell wall mm and bundling with a preferred 
curvature [35j . We assume the latter and write the energy 
for maintaining the crescent shape as 

Ecres = Y J Q du ~ 7r) ’ ( 4 ) 


where u is the arc-length parameter along the crescentin 
bundle attached to the cell wall, c(u) is the local curva¬ 
ture, R c is the preferred radius of curvature, l c is the con¬ 
tour length, and k c is the linear bending rigidity. Equa¬ 
tion (S. 22 ) accounts for the compressive stresses gener¬ 
ated by the crescentin bundle on one side of the cell wall, 
leading to a reduced rate of cell growth, according to 
equation ([I]). As a result, the cell wall grows differen¬ 
tially and maintains a non-zero curvature of the center- 
line. In the absence of crescentin ( k c = 0), our model 
predicts an exponential decay in the cell curvature that 
leads to a straight morphology, consistent with previous 
observations m 

Finally, one must also account for the energy driving 
cell wall constriction. Constriction proceeds via inser¬ 
tion of new peptidoglycan material at the constriction 
site. This process leads to the formation of daughter 
pole caps [35]. We take constriction to be governed by 
an energy of the form E an, = —A S, where S is the surface 
area of the septal cell wall, and A is the energy per unit 
area released during peptidoglycan insertion. 

There exists an optimal cell geometry for a 
given mechanical energy. To apply the model in¬ 
troduced above (equations ([l]) and <§) to interpreting 
the data in Fig. 2, we assume a minimal cell geome¬ 
try given by a toroidal segment with uniform radius of 
curvature R(4>), uniform cross-sectional width w(<j>) and 
the spanning angle To this end, we estimate as 
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FIG. 2. Dynamics of cell shape parameters, (a) Length of the cell medial axis (data shown in black and exponential fit 
from our mechanical model in red), (b) Radius of curvature of the cell medial axis, obtained by calculating the best-fit circle 
to the entire cell (black points), with a time-averaged mean (R) = 4.44±2.12 /im (0 < <j> < 0.5). The mean-field model predicts 
a constant steady-state value for {R((f>)) (red solid line), whereas by accounting for constriction dynamics, the model captures 
the dip in {R(<j>)) seen for 0.5 < <j> < 0.9 (blue solid line), (c) Characteristic cell width, obtained by spatially averaging the 
width at each time point (black points), with a time-averaged mean (w) = 0.74±0.02 /im (0 < <j> < 0.5). The mean field model 
predicts a constant steady-state value for ( w((j >)) (red solid line), whereas cell constriction accounts for the dip in (w(<j))} seen 
for <j> > 0.5 (blue solid line). The shaded regions represent ±1 standard deviation. Model parameters: P = 0.3 MPa, 7 = 50 
nN/pm, k c = 2 nN/im 2 , R c = 0.5 /im, k m = 40 nNgm, R m = 0.31 /im, r = 73 min. 


many mechanical parameters as we can from the liter¬ 
ature and then determine the rest by fitting our exper¬ 
imentally measured values. Turgor pressure in Gram- 
negative bacteria has been measured to be in the range 
0.03 — 0.5 MPa pTTTfTTiTI . We use a value for the effec¬ 
tive internal pressure close to the higher end of the mea¬ 
sured values for turgor pressure, P = 0.3 MPa, in or¬ 
der to account for peptidoglycan insertion. We estimate 
the surface tension as 7 = 50 nN//jm (see Supplemen¬ 
tary model section) and multiply it by the cell surface 
area A(cf>) = nwRd to obtain the cell wall surface energy. 
First, we neglect cell constriction (setting E^iv = 0) and 
assume that the crescentin structure spans the length of 
the cell wall (excluding the endcaps) [40) , with a contour 
length l c {<f>) = (R — w/2)9. The mechanical properties 
of MreB and crescentin are likely similar to those of F- 
actin and intermediate filaments, respectively mm- 
However, due to a lack of direct measurements, we ob¬ 
tain the mechanical parameters k m and k c by fitting the 
model to the experimental data. 

As desired, we find that the total energy E has a sta¬ 
ble absolute minimum at particular values of the cross- 
section diameter w and the centerline radius R, given by 
solution of dE/dw = dE/dR = 0 (see Supplementary 
Fig. 4). The measured values are (w) = 0.74 ± 0.02 /im 
and (R) = 4.44±2.12 /tm (0 < 0 < 0.5), and, as indicated 
by the red solid curves in Fig. 2b,c, the model reproduces 
them with k m = 40 nN/im and k c = 2 nN/nn 2 . While 
the fitted value for k c is numerically close to the estimate 
based on the known mechanical properties of intermedi¬ 
ate filaments (~ 1.5 nN/im 2 ), the value for k m is much 
higher than the bending rigidity of MreB bundles (see 
Supplementary Information). This indicates that k m is 
only determined in part by MreB and can have contribu¬ 
tions from the cell wall. 

Given stable values for w and R , growth is completely 


described by the dynamics of the angle variable 
Consequently, we write the total energy in the scaling 
form E(w, R,0((/))) = 0{<j>)U{w,R), with U the energy 
density along the longitudinal direction. The condition 
for growth then becomes U < 0, such that the energy 
is minimized for increasing values of $(</>). From our ex¬ 
perimental data, the angle spanned by the cell centerline 
increases by an amount A 6 — 0.49 during the entire cy¬ 
cle. Using our parameter estimates and fitting the data 
in Fig. 2, we obtain a numerical value for the energy den¬ 
sity U ~ — 5 nN/tm. We relate the angle dynamics to the 
length by 

J1 = KTl (<t>) = (5) 

acp 7 tw 

where k = — (U < 0) is the rate of longitudinal 
growth, which can be interpreted as resulting from re¬ 
modeling of peptidoglycan subunits with a mean current 
k/ttw, across the cell surface area From an expo¬ 

nential fit to the data for cell length (Fig. 2a), we ob¬ 
tain = 1.6 x 10 - ' 1 (nN /tm min) -1 , which gives us 
an estimate of the friction coefficient, 1/(tv iuR&g) ~ 61 
nN /tm -1 min, associated with longitudinal growth; e.g., 
MreB motion that is known to correlate strongly with 
the insertion of peptidoglycan strands [8] . Our results 
are consistent with previous observations of C. crescen- 
tus cells with arrested division but continued growth 1421 . 

Constriction begins early and proceeds with the 
same time constant as exponential growth. Hav¬ 
ing characterized the dynamics of growth, we now turn to 
constriction at the division plane. As mentioned above, 
we obtain the experimental width at each point along 
each cell’s medial axis. The typical width profile is non- 
uniform along its length, exhibiting a pronounced invagi¬ 
nation near the cell center (with width w m i„ (</>); Fig. lc). 
This invagination, which ultimately becomes the division 
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FIG. 3. Timing and location of th e div ision plane, (a) Time-dependence of (wmin((p)), with the experimental points 
in black and the model fit (equations (S. 131 and |t]|) in red. Fit values: w ma x = 0.805 /im, k^ 1 = 130.92 min, k o = 0.016 
/rm 2 /min. Inset: Minimal geometry of a constricting cell, where S (blue) is the septal cell wall synthesized during constriction, 8 
is longitudinal width of the constriction zone, and ip is the tangent angle at constriction, (b) Ratios of the length of the stalked 
(red) and swarmer (blue) portions divided by the total length (( l st (p)/l((p )) and (l aw (</))/1(<p)), respectively), are approximately 
constant over the cell cycle, (c) Experimental width profile of C. crescentus cells in the initial stage of the cycle ( </> = 0) 
after ensemble-averaging over all data. The width ( w((f> = 0,u)) is plotted as a function of the distance from the stalked end 
normalized to the length of the cell, u/l. Inset: A representative single cell contour immediately after division (p = 0), showing 
the location of the local minimum in the width (wmin, purple line) as well as two local maxima (w^ ax , red line and blue 

line). These two local maxima in the width define the pole regions (shaded in red and blue, respectively), (d) Model width 
profiles of the cell showing symmetric and asymmetric location of the invagination near the midplane for different values of the 
ratio 7 pVt p W =1 (blue, dashed), 1.05 (red, solid) and 1.09 (green, dotted). Parameter values are the same as in Fig. 2 with 
7 p‘ = 5 7 m- The shaded regions in a, b and c represent ±1 standard deviation. 


plane, is readily identifiable early in the cell cycle, even 
before noticeable constriction occurs. We discuss the ki¬ 
netics of constriction in this section, and focus on its loca¬ 
tion later in the manuscript. As shown in Fig. 3a (black 
points), (ui m i n (</>)) progressively decreases towards zero 
until pinching off at cf> = 1. Due to the limited spatial 
resolution of our imaging (phase contrast microscopy), 
the pinch-off process occurring for cp > 0.9 could not be 
captured, but (w m i n (</ > )) at earlier times (i.e., <p < 0.9) is 
precisely determined as a function of </>. 


To model the dynamics of constriction, we assume as in 
Ref. (Fig- 3a, inset): (i) the shape of the zone of con¬ 
striction is given by two intersecting and partially formed 
hemispheres with radii w m ax/2; and (ii) constriction pro¬ 
ceeds by completing the missing parts of the hemisphere 
such that the newly formed cell wall surface maintains 
the curvature of the pre-formed spherical segments. As 
a result, a simple geometric formula is obtained that re¬ 
lates the width of the constriction zone, w m in(</>), to the 


surface area S((f>) of the newly formed cell wall, 


nin (d) — ^maxyl (.S'fc)) /.S) 


max) ) 


( 6 ) 


where S max = i rw^ ax is the maximum surface area 
achieved by the caps as the constriction process is com¬ 
pleted, i.e., when w m i n ((p = 1) = 0. We assume that the 
addition of new cell wall near the division plane initiates 
with a rate, kq, and thereafter grows exponentially with 
a rate, Kd, according to, 

- < ^r = K d S(<p) + K o, (7) 

r dip 


subject to the initial condition S(cp = 0) =0. The first 
term on the right-hand side of equation ([7| follows from 
equation ([lj, using S(<p) as the shape variable, after in¬ 
corporating the constriction energy Edi v (<p). The rate of 
septal peptidoglycan synthesis, Kd, is thus directly pro¬ 
portional to the energy per unit area released during con¬ 
striction, A. The solution, S(4>) = Ko(e KdT,l> — 1 )/iid, 
can then be substituted into equation (S. 13) to derive 
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the time-dependence of w m i n (<j)), whose dynamics is con¬ 
trolled by two time scales: k^ 1 and S max K,Q 1 . 

Fitting equation (S.13) with the data for (w m i n (</>)), we 


obtain nf — 131 min and S max KQ ~ 118 min. The fit¬ 
ted values for the time constants controlling constriction 
dynamics (re^ 1 and S’max^o^ 1 ) are remarkably similar to 
that of exponential cell elongation ((re) -1 — 125 min). 
This shows that septal growth proceeds at a rate com¬ 
parable to longitudinal growth. Therefore, one of the 
main conclusions that we draw is that cell wall constric¬ 
tion (Fig. 3a) is controlled by the same time constant as 
exponential longitudinal growth (Fig. 2a). 

Having determined the dynamics of w m i n (</>), we com¬ 
pute the average width across the entire cell w{4>) using 
the simplified shape of the constriction zone as shown 
in Fig. 3a (inset). The resultant prediction (blue solid 
curve in Fig. 2c) is in excellent agreement with the ex¬ 
perimental data and captures the dip in ( w((f >)) seen for 
cj> > 0.5. Constriction also leads to a drop in the aver¬ 
age radius of curvature of the centerline, as shown by the 
experimental data in Fig. 2b. In the supplementary ma¬ 
terial we derive a relation between the centerline radius of 
curvature R(4>) and the minimum width w m i n (4>), given 
by R~ 1 (dR/d(j>) = ir - in (dtr m i n / d(f>), predicting that cell 
curvature increases at the same rate as w m i n (<j>) drops. 
Using this relation, we are able to quantitatively capture 
the dip in ( R{cf >)) seen for <f> > 0.5 (solid blue curve in Fig. 
2b) without invoking any additional fitting parameters. 

Origin of the asymmetric location of the pri¬ 
mary invagination. We now consider the position of 
the division plane and its interplay with cell shape. As 
shown in Fig. 3b, the distance of the width minimum 
from the stalked pole (I st (<f>)) increases through the cell 
cycle at the same rate as the full length of the growing cell 
(£(</>)), such that their ratio remains constant with time- 
averaged mean (I st /l) = 0.54 ±0.05. The presence of the 
primary invagination early in the cell cycle is reiterated 
in Fig. 3c, which shows the width profile constructed by 
ensemble-averaging over each cell at the timepoint imme¬ 
diately following division. In addition to the width mini¬ 
mum ui m i n (</>), there are two characteristic maxima near 
either pole, w^ ax (</>) and respectively (Fig. 3c, 

inset). As evident in Fig. 3c, the stalked pole diameter 
( u 'max( ( ^ 1 )) is on average larger than its swarmer counter¬ 
part (w^ x (</f>)) (also see Supplementary Fig. 2a). 

We show that the asymmetric location of the invagi¬ 
nation (and the asymmetric width profile) can originate 
from the distinct mechanical properties inherent to the 
pole caps in C. crescentus. The shapes of the cell poles 
can be explained by Laplace’s law that relates the pres¬ 
sure difference, P, across the cell wall to the surface ten¬ 
sions in the stalked or the swarmer pole, r )^ sw . The radii 
of curvature of the poles then follow from Laplace’s law 


2l st, sw 


JDSt.SW 

~ p 


( 8 ) 


where the superscript (si, sw) denotes the stalked or 
the swarmer pole. Thus a larger radius of curvature 
in the poles has to be compensated by a higher sur¬ 
face tension to maintain a constant pressure difference 
P. Assuming that the poles form hemispheres, we have 
P * 4,sw = wfff™ / 2. Our data indicate that the early time 
ratio for 'U’max( < /0/ u 'niax( < /’) (</> < 0.1) shows a strong pos¬ 
itive correlation with the ratio I st (<f>) /l sw (<t>), with an av¬ 
erage value (Wmax/ w mLx) ~ 1.04 (see Supplementary Fig 
3a). Laplace’s law then requires that the stalked pole be 
mechanically stiffer than the swarmer pole; y® 4 > y* 1 ". 
This observation suggests that the asymmetry in the 
lengths of the stalked and swarmer parts of the cell de¬ 
pends upon different mechanical properties of the respec¬ 
tive poles. 

To quantitatively support this claim, we investigate an 
effective contour model for the cell shape. To this end, 
we assume that the fluctuations in cell shape relax more 
rapidly than the time scale of growth. This separation 
of timescales allows us to derive the equation govern¬ 
ing the cell contour by minimizing the total mechanical 
energy (equation §). From the solution we compute 
the resultant width profile for the entire cell (see Supple¬ 
mentary model section). As shown in Fig. 3d, the model 
with asymmetric surface tensions of the poles causes the 
primary invagination to occur away from the cell mid¬ 
plane. The spatial location of the invagination relative 
to the cell length depends linearly on the ratio 7 ^/ 7 ^“ • 
Symmetry is restored for y^/y^ 1 " = 1, as shown in Fig. 
3d (blue dashed curve). We note that a gradient in 7 
along the cell body would imply differences in longitu¬ 
dinal growth rates between the stalked and the swarmer 
portions of the cell (Eq. 0). Our data exclude this pos¬ 
sibility since both I st (4>) and l sw ((f>) grow at the same 
rate re, as evidenced by the constancy of their ratio (Fig. 
3b and Supplementary Fig. 2d). Because C. crescentus 
does not exhibit polar growth, the polar stiffness model 
is consistent with the observed uniformity in longitudi¬ 
nal growth rate. In addition, the non-uniformity in cell 
width comes from the differences in mechanical response 
in the cell wall due to preferential attachment of cres- 
centin along the concave sidewall. For a creS mutant 
cell (where k c = 0), our model predicts a uniform width 
profile before the onset of constriction. 

Cell shape evolution during wall constriction. 
The experimental width profiles show that the growing 
and constricting cells typically develop a second mini¬ 
mum in width (Fig. 4a,b). These secondary invaginations 
are observed in both the stalked and swarmer portions of 
single cells in the predivisional stage (</> > 0.6), although 
they are more common in the stalked portions (Fig. 5a). 
We show here that these secondary minima become the 
primary minima in each of the daughter cells. To study 
the dynamics of the development of the secondary mini¬ 
mum we introduce a new quantity, l^ in (<j>), defined as the 
distance from the stalked pole to the secondary minimum 
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FIG. 4. Comparison of experimental and model cell contours and width profiles, (a) Splined contours of a growing 
and constricting cell at different values of normalized time (j> = 0.0 (red), 0.33 (orange), 0.67 (green) and 1.0 (blue), (b) 
Experimental width profiles plotted against absolute distance from the stalked pole, corresponding to contours in a. (c) 
Contours computed from the cell shape model at different values of w m m and l corresponding to the time points in a. (d) 
Model width profiles corresponding to the contours in c. 


in the stalked part (see Fig. 5c, inset). We find that the 
ratio (^min( < ^)/^ st ( < / > )) h &s a mean value of ~ 0.55 at later 
points in the cell cycle (Fig. 5b), equal to the constant 
ratio maintained by the distance from the stalked pole 
to the primary minimum, (I st (</>))• In fact the ky¬ 
mograph of width profiles (shown over 2 generations for 
a representative single cell) in Fig. 5c demonstrates that 
the predivisional secondary invaginations are inherited 
as primary invaginations after division. This mechanism 
provides continuity and inheritance of the invaginations 
across generations and is an intrinsic element of the mech¬ 
anism for cell division in C. crescentus. 


To quantitatively explain the experimental width pro¬ 
files during constriction, we use our mechanical model 
to determine the instantaneous cell shape by minimiz¬ 
ing the total energy (equation §) at the specified time 
points (see Supplementary model section). To take con¬ 
striction into account, we impose the constraint that 
w(l st ,(j)) = w m i n ((f>), where w m i n (4>) is determined by 
equations (S.13) and (J7|. In addition, we assume non- 
uniform materials properties in the cell wall by taking 
the tension in the cell poles ("fp t ’ sw ) and the septal re¬ 
gion to be higher than the rest of the cell. As constric¬ 
tion proceeds and w m in(</ > ) decreases, we compute the 
shape of the cell contours (Fig. 4c) and the correspond¬ 
ing width profiles (Fig. 4d). The computed width profiles 
faithfully reproduce the secondary invaginations, which 
become more pronounced as the daughter pole caps be¬ 
come prominent. An example of the experimental width 
profiles is shown in Fig. 4b at evenly-spaced intervals in 


time for a single generation, and the corresponding model 
width profiles are shown in Fig. 4d. 

We note that the experimental cell contours in the pre¬ 
divisional stage (cj) > 0.9) bend away from the initial 
midline axis and develop an alternate growth direction 
(Fig. 4a, blue contour). These bend deformations are in¬ 
duced by the microfluidic flow about the pinch-off plane; 
the cells become increasingly “floppy” as the constriction 
proceeds. 


DISCUSSION AND CONCLUSIONS 

The consistent propagation of a specific shape through 
the processes of growth and division relies upon an in¬ 
tricate interplay between the controlled spatiotemporal 
expression and localization of proteins, and cytoskeletal 
structural elements. The high statistical precision of our 
measurements allows us to gain new insights into cell 
morphology. From precise determination of cell contours 
over time, we observe that a typical cell width profile 
is non-uniform at all times with a pronounced primary 
invagination appearing during the earliest stages of the 
cell cycle. During cell constriction, the decrease in the 
minimum width is governed by the same time constant 
as exponential axial growth (Fig. 3a). Furthermore, the 
location of the primary invagination divides the cell con¬ 
tour into its stalked and swarmer compartments, such 
that the ratio of the length of the stalked part I st (<f>) to 
the total pole-to-pole length l(4>) remains constant dur- 
















FIG. 5. Location of division plane is set in the previous generation, (a) Value of the minimum width normalized by 
the respective maximum width for the stalked ({Wnin(^)/tBmaxW)) and the swarmer (((/>)/w^, x (0))) parts, (b) Ratio of 
length from stalked pole to secondary minimum normalized by length from stalked pole to primary minimum / I st (<p)), 

green) and ratio of length from stalked pole to primary minimum normalized by total length {{l st {<t>)/l{<t>)), red). The former 
ratio remains constant at ~ 0.55, while the latter obtains this value at the end of the cell cycle. In comparing averages between 
two generations here, we indicate values of <j> from the first generation as negative (e.g., f> = —0.2 is 20% of the way from the 
subsequent division). The green points are not shown for <j> < —0.25 due to increased errors in identification, (c) Kymograph 
of width profiles for a typical cell over two generations. The time evolution of the widths (color scale) illustrates continuity of 
the location of the minima across generations. That is, the location of the secondary minimum just before division (Z^m, white 
line) becomes the primary minimum ( I st , black line) just after division (horizontal dashed line). The schematic at right shows 
two measured contours that correspond to time points immediately before = —0.05) and after (cp = +0.05) the division 
event shown in the kymograph. 


ing the cycle with a mean value {I st (cp) / l((p)) ~ 0.55 (Fig. 
3b). These observations and our mechanical model lead 
to two important conclusions: first, the dynamics of cell 
wall constriction and septal growth occur concomitantly , 
and second, the asymmetric location of the primary in¬ 
vagination can be explained by the differences in mechan¬ 
ical properties in the stalked and swarmer poles. A corol¬ 
lary of the first conclusion is that the size ratio threshold 
at division occurs naturally without requiring a complex 
timing mechanism m- 

In addition to the primary septal invagination, the 
cell contours exhibit a pronounced secondary invagina¬ 
tion during the predivisional stages (Fig. 4). Remark¬ 
ably, the secondary invaginations develop at a precise 
location relative to the total length of the stalked com¬ 
partments, {l^ in (cp)/I st ((/))} ~ 0.55 (Fig. 5b). The data 
thus allow a third conclusion: these secondary invagina¬ 
tions are inherited as primary invaginations in each of 
the daughter cells, directing the formation of the division 
plane in the next generation. Thus, through consistent 
and controlled nucleation of invaginations across genera¬ 
tions, C. crescentus cells maintain a constant ratio of the 


sizes of stalked and swarmer daughter cells. 

Our experimental observations and the parameters in 
the cell shape model can be related to the current molec¬ 
ular understanding for Gram-negative bacteria, in par¬ 
ticular C. crescentus. Before the onset of noticeable 
constriction, cell shape is dictated by the mechanical 
properties of the peptidoglycan cell wall in addition to 
various shape-controlling proteins such as MreB, MreC, 
RodZ and CreS. Single molecule tracking studies have 
revealed that MreB forms short filamentous bundles an¬ 
chored to the inner surface of the cell wall and moves 
circumferentially at a rate much faster than the rate of 
cell growth Emu- In vitro experiments show that MreB 
filaments can induce indentation of lipid membranes, sug¬ 
gesting that they may have a preferred radius of curva¬ 
ture [46!. Thus on time scales comparable to cell growth, 
.Ewidth is determined in part by the energy cost of ad¬ 
hering MreB bundles to the cell wall (see Supplementary 
model section). 

Bacterial cell division is driven by a large complex of 
proteins, commonly known as divisomes that assemble 
into the Z-ring structure near the longitudinal mid-plane 
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of the cell m The Z-ring contains FtsZ protofilaments 
that are assembled in a patchy band-like structure [47] . 
FtsZ protofilaments are anchored to the cell membrane 
via FtsA and ZipA, and play a crucial role in driving 
cell wall constriction [3B|- During constriction, the divi- 
some proteins also control peptidoglycan synthesis and 
direct the formation of new cell wall via the activity of 
penicillin-binding proteins (PBPs) [49] [5(J. Thus the di- 
visome plays a two-fold role by concomitantly guiding 
cell wall constriction and growth of the septal peptido¬ 
glycan layer. According to our model the constriction of 
the cell wall is driven by the synthesis of septal cell wall 
at a rate nd (~ (k)), which can be directly related to 
the activity of PBPs triggered by the divisome assembly. 
Furthermore, in our model it is sufficient that the divi¬ 
some guide the curvature of cell wall growth in the septal 
region (see Fig. 3a, inset). 

While the mechanism behind the precise asymmetric 
location of the division plane in C. crescentus cells is not 
well understood, it is likely that the ATPase MipZ helps 
division site placement by exhibiting an asymmetric con¬ 
centration gradient during the predivisional stage m- 
MipZ activity inhibits FtsZ assembly; as a result of polar 
localization of MipZ, Z-ring assembly is promoted near 
the mid-cell [52]. Our cell shape model suggests that the 
early time asymmetric location of the primary invagina¬ 
tion, which develops into the division plane, is controlled 
by the differences in surface tensions maintained in the 
poles. The presence of this invagination at <f> = 0, as 
inherited from the secondary invaginations in the pre¬ 
vious generation, aids in Z-ring assembly at the site of 
the invagination. The curvature-sensing capability of the 
Z-ring may be enabled by the minimization of the FtsZ 
polymer conformational energy that is determined by the 
difference between cell surface curvature and FtsZ spon¬ 
taneous curvature mm- 

A higher tension in the stalked pole can be induced 
by asymmetric localization of polar proteins, such as 
PopZ, early in the cell cycle. Experiments have shown 
that PopZ localizes to the stalked pole during the initial 
phase of the cell cycle and increasingly accumulates at the 
swarmer pole as the cell cycle proceeds [Mj. Consistent 
with this observation, our data show that the correlation 
between the pole sizes (determined by the ratio of sur¬ 
face tension to pressure) and the stalked and swarmer 
compartment lengths tend to disappear later in the cycle 
(Supplementary Fig. 3), as cell constriction proceeds. A 
recent experimental study also demonstrates that molec¬ 
ular perturbation of Clp proteases can destroy the asym¬ 
metry of cell division in C. crescentus [55], suggesting the 
interplay of subcellular protease activity with the physi¬ 
cal properties of the cell wall. 

Earlier theoretical models have predicted that a small 
amount of pincli-off force from the Z-ring (~ 8 pN) is 
sufficient to accomplish division by establishing a direc¬ 
tion along which new peptidoglycan strands can be in¬ 


serted [271 . In contrast, our data combined with the 
mathematical model allows the interpretation that the 
early time asymmetric invagination in the cell wall can 
set the direction for the insertion of new peptidoglycan 
strands. Constriction results from exponential growth of 
surface area in the septum (at the same rate as longi¬ 
tudinal extension). The instantaneous cell shape is de¬ 
termined by minimizing the energy functional at given 
values of the cell size parameters. 

Finally, from our estimate of the cell wall energy den¬ 
sity U (~ — 5 nN/jm), we predict that a net amount 
A0|E/| ~ 2.5 nN/zm of mechanical energy is used by the 
peptidoglycan network for cell wall growth. For a C. cres¬ 
centus cell of surface area 12.5 —25^zm 2 , layered with gly- 
can strands of length ~5 nm and cross-linked by peptide 
chains with maximally stretched length ~4 nm [5], there 
are roughly 10 6 peptidoglycan subunits. Thus on aver¬ 
age, each peptidoglycan subunit can consume mechanical 
energy of ^2.4xl0 -6 nN/zm, or ~0.6 ksT at a temper¬ 
ature T = 31°C. Cell wall remodeling and insertion of 
new peptidoglycan material can likely create defects in 
the peptidoglycan network [56]. One thus expects cellu¬ 
lar materials properties to change over time, as a result of 
these molecular scale fluctuations. Although we neglect 
such variations in our mean field model, it nonetheless 
quantitatively captures the average trends in cell shape 
features. In future work we plan to more closely connect 
the energy terms of the continuum model with molecular 
details. 


METHODS 

Acquisition of experimental data. Data were ac¬ 
quired as in Ref. mi- Briefly, the inducibly-sticky 
Caulobacter crescentus strain FC1428 was introduced 
into a microfluidic device and cells were incubated for 
one hour in the presence of the vanillate inducer. The 
device was placed inside a homemade acrylic microscope 
enclosure (39" x 28" x 27") equilibrated to 31°C (tem¬ 
perature controller: CSC32J, Omega and heater fan: 
HGL419, Omega). At the start of the experiment, com¬ 
plex medium (peptone-yeast extract; PYE) was infused 
through the channel at a constant flow rate of 7 /zL/min 
(PHD2000, Harvard Apparatus), which flushed out non¬ 
adherent cells. A microscope (Nikon Ti Eclipse with the 
“perfect focus” system) and robotic XY stage (Prior Sci¬ 
entific ProScan III) under computerized control (Lab- 
View 8.6, National Instrument) were used to acquire 
phase-contrast images at a magnification of 250X (EM- 
CCD: Andor iXon+ DU888 lk x lk pixels, objective: 
Nikon Plan Fluor 100X oil objective plus 2.5X expander, 
lamp: Nikon C-HFGI) and a frame rate of 1 frame/min 
for 15 unique fields of view over 48 hours. In the present 
study we use a dataset consisting of 260 cells, correspond¬ 
ing to 9672 generations (division events). 
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Analysis of single cell shape. The acquired phase- 
contrast images were analyzed using a novel routine we 
developed (written in Python). Each image was pro¬ 
cessed with a pixel-based edge detection algorithm that 
applied a local smoothing filter, followed by a bottom-hat 
operation. The boundary of each cell was identified by 
thresholding the filtered image. A smoothing B-spline 
was interpolated through the boundary pixels to con¬ 
struct each cell contour. Each identified cell was then 
tracked over time to build a full time series. We chose to 
include only cells that divided for more than 10 genera¬ 
tions in the analysis. A minimal amount of filtering was 
applied to each growth curve to remove spurious points 
(e.g., resulting from cells coming together and touching, 
or cells twisting out of plane). The timing of every divi¬ 
sion was verified by visual inspection of the corresponding 
phase contrast images, so that the error in this quantity 
is approximately set by the image acquisition rate of 1 
frame/min. 
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Supplemental Material 

EXPERIMENTAL CELL SHAPE PARAMETERS 


A number of quantities can be immediately calculated from each splined cell contour (Supplementary Fig. la), 
including the cross-sectional area. The cell medial axis was determined by calculating the Voronoi diagram of the cell 
contour m , pruning the branches (Supplementary Fig. lb), and extending this skeleton to the edges of the cell contour 
in a manner that preserves average curvature of the medial axis (Supplementary Fig. lc). The intersections between 
the cell medial axis and contour represent the stalked and swarmer poles, respectively. The cell length was calculated 
by evaluating the distance along the medial axis between either pole. The cell widths were determined by creating 
ribs perpendicular to the medial axis along its length and determining the distances between their intersections with 
opposite sides of the cell contour (Supplementary Fig. Id). To calculate time-averaged quantities, we normalized 
trajectories for each generation by the respective division time r (thus converting each variable to a function of 
(f> = t/r). We then split these data into 73 bins, as (r) = 73 ± 7 min under these conditions, and ensemble-averaged 
each of these bins over every generation. 


a c 



012345 012345 

x (pm) x (pm) 


Supplementary Figure 1. Procedure for calculating shape parameters, (a) B-spline cell contour extracted from 
segmented phase contrast image, (b) The interior part of the medial cell axis was obtained by taking the Voronoi diagram of 
the cell contour, which is equivalent to the locus of inscribed circles, (c) The medial axis located on each side of the width 
minimum (excluding the constriction zone) was extended to either pole such that average radius of curvature remained constant, 
(d) The cell widths are the ribs perpendicular to the medial axis. 


We can define five specific values of the width according to the local minima and maxima in the width profile (w m i n , 
w^ax, jjx, Wmin; W min)i n °t °f which may be present in any given cell. Where we identify a primary minimum 
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in the width profile, w m i n , we can also determine two local maxima, ax and corresponding respectively to 

the stalked and swarmer portions (Supplementary Fig. 2a). These values are approximately constant throughout 
the cell cycle, with (u>^ ax (^)) > However, at later times ((/> > 0.6) the value of u^ x (</>) increases 

until ( , u ; max(^)) ~ («Cx('/'))• some cases, additional secondary local minima are observed, w^ in and w^ in , 
corresponding respectively to the stalked and swarmer portions (Supplementary Fig. 2b). Although we note the value 
of these quantities for early times here (where they are approximately equal to their respective local maxima), these 
minima can only be determined with certainty at later times ( (f> > 0.6). There, we observe the presence of a secondary 
minimum in the stalked portions of most cells. 



mean ± S.D. 


</> 

w 

0.74 ± 0.02 

/im 

[0, 0.5] 

Wmax 

0.85 ± 0.02 

/im 

[0,0.5] 

w st - 

lu min 

0.84 ± 0.03 

/im 

[0,0.5] 

w sw 

'"max 

0.83 ± 0.02 

/im 

[0, 0.5] 

yj S ™ 

< "min 

0.82 ± 0.02 

/im 

[0,0.5] 

R 

4.44 ± 2.12 

/im 

[0, 0.5] 

R at 

4.34 ± 2.50 

/im 

[0,0.5] 

RSW 

3.94 ± 2.42 

/im 

[0, 0.5] 

I st /l 

54.3 ± 5.2 

% 

[0,1] 

l sw /l 

45.7 ± 5.2 

% 

[0,1] 


Supplementary Table 1. Average values of the cell shape parameters that are assumed constant over the 
specified time interval of </>. The time-average value of the width ( w ) is smaller than both (w^ in ) and (wj^; n ) because it 
accounts for the width at the primary minimum w m i n . Note the close correspondence between either width maximum and its 
corresponding minimum for (j> < 0.6, as well as u>^ ax and in . Both width and radius of curvature for the swarmer portion are 
smaller than their respective values for the stalked portion. 

We split each cell into stalked and swarmer portions according to the location of the primary minimum in the width 
profile. The radius of curvature ( R ) was calculated as the radius of the best-fit circle to the cell medial axis. We tested 
two methods of determining the radius of curvature: (I) fitting the whole cell to a circle or (II) fitting the stalked and 
swarmer portions separately. At early times in the cell cycle, (II) gives poorer results because of fewer data points. 
At later times in the cell cycle, (I) gives poorer results because the stalked and swarmer portions can indeed have 
differing radii of curvature, which we attribute to the alignment of the swarmer portion of the cell with the direction 
of fluid flow after the division plane has narrowed enough that it becomes mechanically decoupled from the stalked 
portion, i.e., like a flexible hinge. However, in the mean the results of (II) are equal to the value calculated by (I) 
for earlier times. Therefore, we use only data from method (I) but exclude later time points (Supplementary Fig. 
2c). Because of the variation in the values of the calculated radius of curvature was large, such that we cannot define 
a reasonable arithmetic mean without arbitrarily filtering the dataset, we first averaged the corresponding unsigned 
curvature (equivalent to ( R ~ 1 )) and then converted to radius of curvature (i.e., what we report as ( R) is actually the 
harmonic mean, calculated as (I? -1 ) -1 ). 

The length was also split according to the locations of the local minima, into I st (distance along cell medial axis 
from stalked pole to u> m i n ), l sw (distance along cell medial axis from swarmer pole to ic m i n ), Z^ in (distance along cell 
medial axis from stalked pole to u>*( in ), and l^ n (distance along cell medial axis from swarmer pole to u^ n ). The 
value of I st and /) r ( in are compared in Supplementary Fig. 2d. Note that the length of the stalked portion I st grows 
exponentially, with the same time constant as the length l (as does l sw , although it is not shown here for clarity), 
which is a necessary condition for the addition of peptidoglycan material along the entire length of the cell when the 
location of is set at early times. At later times ( (f> > 0.6), the length from stalked pole to secondary minimum 
l^ in also starts to increase. 

We relate the asymmetry in the length of stalked and swarmer portions at early times to asymmetries in the stalked 
and swarmer poles. The model predicts a linear relationship between the ratio of lengths I st /l sw and the ratio of the 
tensions at either pole We cannot directly measure the latter quantity, but from Laplace’s law it is equal 

to the inverse ratio of the mean curvatures at either pole, which we approximate as the inverse of half the maximum 
pole width, assuming that the poles are hemispheres with diameter and ic^!( x , respectively. Supplementary Fig. 
3 shows scatter plots comparing «i„ax/*max to I st /l sw for three different time intervals. The red best-fit line is shown 
only for Supplementary Fig. 3a, for which R 2 = 0.15 (linear fits to all other plots produced values of R 2 < 0.05). This 
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Supplementary Figure 2. Unnormalized values of the cell shape parameters. Time-dependence of (a) width maxima, 
and (b) width minima, of both stalked and swarmer portions. These data correspond to Fig. 5a in the main text, (c) Radius 
of curvature calculated for either the entire cell (black), the stalked portion only (red), or the swarmer portion only (blue), (d) 
Unnormalized length from stalked pole to either the primary (red) or secondary (green) width minima. These data correspond 
to Fig. 5b in the main text. Superscript * refers to either st or sw, respectively, as indicated in each figure legend. 


line runs from the point (1.0,1.01), corresponding to the dashed blue curve in Fig. 3c of the main text (which is the 
case of a symmetric cell) to (1.2,1.04), corresponding to the red solid curve in Fig. 3d of the main text (which is the 
case of the average asymmetric C. crescentus cell). Note that at times <j> > 0.6, any correlation disappears. 

In order to quantify the error in our width profiles, we imaged a single field of view of 24 C. crescentus cells perfused 
in complex medium at 31° C at a frame rate of 5 frames per second (300 times faster than the frame rate used to 
acquire all other data), and calculated the splined contours for each cell. We focused in particular on a single “dead” 
(non-growing and non-dividing) cell, and found the root-mean-square deviation (RMSD) of nearest points along the 
cell contour between subsequent frames to be 12 mn. In addition, we found the RMSD of equivalent points along the 
width profile between subsequent frames to be 28 nm, or a 3.2% pinch depth at the average value of w^ ax = 0.85 ^m. 


CELL SHAPE MODEL 

Cell Wall Mechanics. The total energy E for the bacterial cell wall is given as the sum of contributions from 
an active internal pressure P driving cell volume (V) expansion, mechanical energy E wa ,n in the cell wall, and the 
mechanical energy of interactions with cytoskeletal proteins -Eproteins: 


E — —PV + Inwall + ^proteins • (S.l) 

The bacterial cell wall consists of a network of glycan strands cross-linked by peptide chains know as the peptidoglycan 
network. Growth occurs via the insertion of new peptidoglycan strands into the existing network along with the 
breaking of existing bonds due to turgor pressure induced stretching. We assume that elastic equilibrium is reached 
rapidly as compared to the rate of synthesis of new material m ■ As a result of cell wall remodeling and irreversible 
elongation, growth can be understood as resulting from plastic deformations [aMj- To understand the origin of the 
cell wall tension, 7 , in the model introduced in the main text, we consider the cell wall as a thin elastic shell that 
deforms plastically when stretched beyond a maximum strain £y, the yield strain. A thin shell has two modes of 
elastic deformations, bending and stretching m, such that Inwall = -E-stretch + Abend- In the limit of small thickness 
of the shell h as compared to its radii of curvature, one can neglect the bending energy (that scales as Abend ~ h 3 ) 
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Supplementary Figure 3. Ratios of maximum width at stalked:swarmer poles versus stalked:swarmer lengths. 

The three columns delineate the time period from which data were taken: (a) first 10% of the generation, (b) middle 10% of 
the generation, and (c) last 10% of the generation. Blue points show the average value of irix/Wmax calculated by binning 
I st /l 3W into bins of width 0.05, with size of each error bar showing the standard deviation of the values within that bin. The 
best-fit line for 1 < I st /l sw < 1.2 is indicated in red (coefficient of determination R 2 = 0.15). This line runs from (1.0,1.01) to 
(1.2,1.04). 


whereas the stretching energy E s te t c h is given by, 

J /> 

-^stretch = ^ I (S.2) 

where cr^- is the mechanical stress tensor and is the strain tensor. As yield strain is reached at the onset of growth, 
we have £ij = £ydij (assuming isotropic stretching), where Sij is the Kronecker delta. Furthermore, assuming a 
Hookean constitutive relation for the stress tensor m 

Y 

&ij — ^ 2 [("i T" L / £kk&ij\ 7 (S- 3 ) 


where Y is the Young’s modulus and v is the Poisson ratio, we have F4tretch = jA, where, 

1 Y key 

7 = 2 = ■ 


(S.4) 


For a Gram-negative bacterial cell wall of thickness h ~ 3 nm, elastic modulus Y ~ 40 MPa m and average yield 
strain Sy ~ 0.5, the wall tension is estimated to be 7 = 50 nN//xm. While the actual value for the turgor pressure 
counteracting this tension can be contested, our choice for the numerical value for internal pressure P can be justified 
using a simple mechanical argument. Radial force-balance dictates that in order to maintain an average cross-sectional 
radius r ~ 0.4 /im, a cell wall with surface tension 7 ~ 50 nm//mi has to balance an internal pressure of magnitude 
P = 27 /r — 0.25 MPa, which is numerically very close to our choice for P = 0.3 MPa. 

Cytoskeletal bundle mechanics. Next, we model the mechanical energy in the cell wall due to interactions with 
cytoskeletal proteins. Our minimalist approach considers two crucial protein bundles that are directly responsible 
for maintaining the shape of C. crescentus cells. MreB protein bundles contribute an effective energy i^width, which 
favors a rod-like shape, and crescentin filament bundles contribute an energy E cres , which favors a crescent-like shape. 
We thus have E prote i ns = E w idth + E cres . MreB subunits form patchy filamentous bundles adherent to the cell wall 
and oriented perpendicular to the long axis of the cell [331121. The elastic energy stored in an adherent MreB subunit 
is given by 


AfvireBy 




(S.5) 


where i labels the subunit, r t is the circumferential radius of curvature of the cell wall where the subunit is attached, li 
is the length of the subunit, R m is the intrinsic radius of curvature and kb is the bending rigidity of the associated MreB 
bundle. The total energy imparted by a collection of N rn attached MreB subunits is given by if w idtii = -E’MreBy- 

Next, employing a continuum mean field assumption, we replace individual subunit lengths by their average length 
(£) ~ 5 nm [ 310| and assume a uniform number density p m of MreB subunits in the cell surface to obtain 


— dA - - 


P,i 


Awidth — 


2 


(S.6) 
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where, k m = kbp m (£) is the effective bending modulus due to MreB induced traction forces. The bending rigidity 
of an MreB bundle is given by kb = fcMreB?^, where /cMreB is the flexural rigidity of MreB filaments (assumed to 
be similar to F-actin), n is the number of MreB protofilaments per bundle, and £ is an exponent in the range 1-2 
depending on the strength of crosslinking or bundling agents ■ MreB filaments in a bundle appear to have strong 
lateral interaction with negligible filament sliding, so we assume £ ~ 2. The diameter of an MreB protofilament is 
~ 4 nm im and the average width of an MreB bundle has been determined from super-resolution imaging to be in 
the range 60-90 nm giving the estimate n ~ 15-22. Estimating the surface number density of MreB subunits as 
Pm ~ 5 x 10 4 /im ~ 2 and using fc mre b — 10 ~ 4 nN/jm 2 [3H, we obtain k m in the range 5.6-12.1 nN/irn. 

Crescentin proteins form a cohesive bundled structure anchored to the sidewall of C. crescentus cells. The energetic 
contribution due to crescentin is given by 


E — 

-L^c. res — 


kc 

T 


ds 




(S.7) 


where k c is the bending rigidity, s is an arc-length parameter, c(s) is the longitudinal curvature of the cell wall and 
R c is the preferred radius of curvature of crescentin bundles. The bending rigidity of crescentin can be expressed as 
k c = Y C I , where Y c is the Young’s modulus and I is the area moment of inertia of the bundle (with width w c ) given 
by I = 7tw 4 /64. Since crescentin is an intermediate filament homologue, we assume that Y c is similar to the Young’s 
modulus of intermediate filament bundles given by Y c ~ 300 MPa [323 • Assuming w c ~ 0.1 /im, we estimate the 
bending rigidity of a crescentin bundle to be k c ~ 1.5 nNgm 2 . 

Mean field model for cell shape and size dynamics. As described in the main text, the total mechanical 
energy in the cell wall of C. crescentus cells can be given as a sum of contributions from internal pressure, wall surface 
tension, mechanical energy of interactions with bundles of cytoskeletal proteins such as crescentin and MreB, and the 
constriction energy E& j v during cell division. In the mean field description, we neglect the contributions from E d i v 
and disregard any spatial variations in cell geometry. The mean field description is a good approximation of the cell 
shape dynamics for <fi < 0.5, where the average width and the midline radius of curvature remains constant (Fig. 2 in 
the main text). We approximate the shape of C. crescentus cells as the segment of a torus with radius of curvature 
R{t), cross-sectional width w(t) and spanning angle 9{t). We also neglect the pole caps that are mechanically rigid 
and do not remodel during wall growth. The total energy is then given by 

E[w, R, 9] = -P(ttw 2 R9)/4 + y(7r RwO) + £ w idth + Acres , (S.8) 


where 

A'width = —~- (ttRw9) [{w/2)~ l - R- 1 } 2 , 


(S.9) 


Acres = y (R ~ w/2)0 [(I? - w/ 2)" 1 - R~ l f . 


(S.10) 


From the above expressions, we see that the total energy has the scaling form E[w,R,6\ = 9U[w,R]. According to 
Eq. (1) in the main text, the dynamics of the radius of curvature R, the spanning angle 9 and the width w are given 
by 


1 dR _ dE 
R~dt ~ ~ <i>R dR ’ 
1 dw dE 

w dt w dui 

1 d9 dE 

~9~dt ~ ^ e ~d9 ’ 


(S.lla) 

(S.llb) 

(S.llc) 


where and are the rate constants. The steady-state values for the radius of curvature of the centerline 

and the cross-sectional width is given by the solutions to dnU = 0 and d w U = 0. As shown in Supplementary Fig. 
4a the energy density U admits a stable absolute minimum in width, controlled by the stiffness parameter k m such 
that for lower values of k m (red curve in Supplementary Fig. 4a), U does not have a minimum and the width grows in 
time rendering the rod-like shape unstable. Similarly, the parameter k c controls cell curvature, such that the energy 
density has a stable absolute minimum in R beyond a critical stiffness k c (Supplementary Fig. 4b). At lower values 
of k c the cell does not have a stable radius of curvature and the energy is minimized by increasing R. This leads to 
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straightening of the cell’s medial axis. In the parameter range where the cell maintains a stable value for w and R, 
the energy density U becomes a constant during cell elongation. From the growth law introduced in the main text, 
the spanning angle 9 evolves in time according to the equation 


f = -c t.um , 


such that the condition for exponential growth becomes U < 0. 



w (jum) 



(S.12) 


Supplementary Figure 4. Dependence of the mechanical energy on cell width and radius of curvature, (a) 

Energy density U as a function of mean cell width w at various values of bending stiffness (in units of nlVpm): k m = 30 (red), 
k m = 40 (green) and km = 50 (blue), for a fixed radius of curvature R = 4.4 /im and k c = 2 nN/im 2 . The absolute minima (if 
it exists) are indicated by solid circles, (b) U as a function of radius of curvature R at various values of crescentin stiffness (in 
units of nN (im 2 ): k c = 1.9 (red), k c = 2.0 (green) and k c = 2.1 (blue), for a fixed width w = 0.74 /im and km = 40 nN/tm. 
Other parameters: P = 0.3 MPa, 7 = 50 nN/ //m, R c = 0.5 /tm, R m = 0.31 /mi. 


Constriction dynamics. As described in the main text, we assume that the shape of the constriction zone is 
given by two intersecting hemishperical segments with diameter w max . The total surface area of the septum is given 
by S(t) = 7 Tw max l B (t), where l s (t) is the total length of the spherical segments (see Supplementary Fig. 5a). Using 
elementary geometry, we obtain the following relation between l s (t) and the minimum width, w m i n (t): 

TCrriin (t) — 7C max \/l {Is (t)/w max )2 . (S.13) 

If septal growth occurs by addition of new peptidoglycan strands while maintaining the curvature of the spherical 
segments, then l s (t) is the shape variable controlling the growth of septal surface area S as well as the constriction 
dynamics of w rn i n (t). To describe the dynamics of l s (t ), we assume that initial phase of elongation occurs with a 
velocity Vo, and thereafter growth follows an exponential law with rate n d , 

-^r= v o + K d l s (t). (S.14) 

Multiplying both sides of the equation by 7 rw max , we derive the growth dynamics of the surface area S(t), 

^ = k 0 + K d S{t) , (S.15) 

where Kq = 7 rw max r! 0 . Having determined the dynamics of w m i n , one can evaluate the time-dependence of the average 
width w defined as, w(t) = w(u,t)du, where u is the coordinate along the centerline. Using the simplified 

geometry of the constricting cell, given by two toroidal segments (peripheral regions) connected by two intersecting 
hemispheres (septal region), one gets, 


)(f) = 


m 


w(l(t) - l a (t)) 


■ M/2 


5 (it, t)dv 


(S.16) 


where w s (u) = w max y/l — ( 2 it/u; max ) 2 , is the width in the septal region and w is the average width of the stalked and 


swarmer components as determined by Eq. (S.llb). As shown in Fig. 2c in the main text (blue solid line), Eq. (S.16) 
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is in excellent agreement with the experimental data and captures the dip in w(t) seen for </> > 0.5. Fig. 2b in the 
main text also shows that the radius of curvature of the centerline ( R(t )) drops for cf> > 0.5 when cell constriction is 
prominent. Therefore we examine the role of constriction dynamics on the time-dependence of the centerline radius 
of curvature. 

The centerline spans a distance l s (t ) and an angle 0 s {t) in the septal region, as shown in Supplementary Fig. 5a. 
Consider a small segment 5l s {t ) along the centerline with a local radius of curvature R and spanning angle 50 s . We 
then have the relation 5l s (t) = R(t)56 s (t). This leads to the following identity, 


1 d(5l s ) 1 dR 1 d{50 s ) 

51 s dt R dt + 50 s dt 


Now the geometry of the septal region (Supplementary Fig. 5a) directly relates the rate of increase in the spanning angle 
50 s {t) to the rate of drop in the tangent angle (ip(t)) at the constriction site, jg- ■ Furthermore, using 

the relations, w max sin ip(t) = w m in(i) and w max cosipit) = l s (t), we get 5ip(t) = 5w min {t)/l s (t) = -5l s (t)/w m i n (t). 


The last equality follows after variations of Eq. (S.13). As a result we have the following kinetic relation. 


d(5ip)/dt dw m i n /dt d(5l s ) d(56 s )/dt 


Sip 


5L 


50s 


Combining the identities in Eq. (S.17) and (S.18) we get, R 1 dR/dt = w m [ n dw, 


(S.18) 

i/dt, showing that the radius of 


curvature drops at the same rate as w m i n (t) shrinks. 




Supplementary Figure 5. Geometry of a constricting cell contour, (a) Schematic of the constriction zone of a dividing 
cell where the dashed line indicates the centerline with length l s , radius of curvature R and spanning angle 0 S . S (blue) denotes 
the septal cell wall, ip is the tangent angle at the constriction site and 5 is the width of the division ring, (b) Compartmentalizing 
the cell contour into the (1) stalked pole p s t (red), (2) swarmer pole p sw (purple) (3) upper wall 1 st / 3W (with crescentin bundle 
attached, green), (4) Lower wall 2 st / sw (black) and the (5) Septal region (blue). 

Contour Model for cell shape. To quantitatively capture the experimentally observed spatial variations in 
cell shape we study an effective two-dimensional contour model for the cell shape. This approach facilitates closer 
comparison with the two-dimensional splined contours obtained from our experimental data. The contour description 
of the rod-like cell shape negelects circumferential variations in cell geometry. The model incorporates non-uniform 
materials properties and mechanical constraints across the cell wall, with the poles and the septal region being 
mechanically stiffer than the rest of the cell wall. By exploiting the rod-like geometry, one can use the centerline 
curve to divide the contour into two parts, the upper and the lower cell wall. As shown in Supplementary Fig. 5b we 
further subdivide the cell contour into the stalked and swarmer poles (p a t/sw)> the upper (1 s t/sw) and lower cell wall 
( 2 St/aw), and the septal region (sep). 

We parametrize the instantaneous shape of the cell contour using the two-dimensional centerline curve R(it, f), 
the length of the centerline l(t) and the width w(u,t), where u is the absolute distance along the centerline from 
the stalked pole such that u £ [0, l(t)]. If n denotes the outward unit normal vector on the centerline, the curves 
defining the upper and lower parts of the cell contour, r± (u,t), are given by the relation, r± = R ± w± n, where w± ( u ) 
represent the perpendicular distances of the top and bottom curves from the centerline. The total cell width is then 
given by w(u, t) = w+(u, t) + W- (u, t). 
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It is convenient to switch to polar coordinates, where the shape of the cell contour is given by the re-parametrized 
curve r ±(<p,t), where 0 < ip < 9(t ) is the angular coordinate spanning the centerline, which can be approximated as 
the arc of a circle with radius R(t). Since the ratio w±/R (~ 0.1) is small at all times, one can approximate the local 
curvature as 


c±(<p,t) 


1 f 1 _ w' ± (<p,t) + w'l(ip,t) 
R(t) \ R(t) 


u ±i UJ ±-> 1 


(S.19) 


where prime denotes derivative with respect to ip and the subscripts ± represent the upper and the lower part of 
the cell contour respectively. Furthermore, in the linear regime, the differential arc length can be approximated as 
du — Rdip. The dynamics of the shape parameters l(t), 0(t), and R(t) are determined from the kinetic law in Eq. (1) 
of the main text. 

The instantaneous width profile w(u,t ) results from minimizing the total energy functional, which leads to the 
following shape equation, 


P - 7 (u)c±(u,t) + /width (u,t) + /cres(^?^) — 0 i (S.20) 

where, 7 (u) is the tension on the cell contour and /width(«, t ) and / cres (+ t) are the linear force densities on the contour 
due to maintenance of width and the crescent shape, respectively. In the contour model we simplify the energetic 
contribution due to maintenance of width (acting on the sections 1 S t/ S w an d 2 s t/sw) as 

T^width = ^ 7 p J du(w± - R m f , (S. 21 ) 

where the elastic constant K m depends on the bending rigidity k m introduced in the main text as K m = k m /Rp n . This 
linear approximation holds if (w± — R m )/Pm +C 1. From our data and mean field model fits we get (w± — R^/Rm ~ 
0.15. The resultant force density is given by /width (u,t) = —K m ( w±(u,t ) — R m )- 

The bending energy induced by crescentin protein bundles anchored onto the cell wall (regions l s t/sw) is given by 

Ecies = y J ds [c + (u) - c 0 ] 2 , (S. 22 ) 

where Co = 1 /R c is the spontaneous curvature of the crescentin bundle, and s is the arc-length parameter along the 
upper part of the cell contour (l s t/ S!U ) which is related to u as ds = du — w + dip. To obtain the force density we 
consider an infinitesimal deformation, <5r, of the upper contour as r + —» r + + n + i5r, where n + is the outward unit 
normal. Accordingly the curvature and the differential arc-length changes, c+ —> C+ + <5c_|_ and ds —> ds + S(ds), where 
Sc + = c 2 + 5r + d 2 (Sr)/ds 2 and S(ds) = —c + ds [fHH], The resultant force density is obtained after variations of the 
energy functional, SE crea = f f cres Srds , where f cies 5rds = £[(c+ — Co) 2 ds]. This leads to the following non-linear force 
contribution: 


/cres 


gFc± 3 2 

ds 2 + + + 0 


Using Eq. (S.191, / cres can be linearized and expressed using the angular coordinate ip as 

1 


fcies{u,t) — 


t) + -( CqR 2 - 5)to+(u, t) + -(cq !? 2 — 3)w+(u, t) + R( 1 - CqI? 2 ) 


(S.23) 


(S.24) 


Compartmentalizing the cell contour. (1) Pole caps: The cell poles are assumed to be mechanically inert in 
the sense that they do not interact with the active cytoskeletal proteins such that /width = /cres = 0. The mechanical 
forces acting on the cell poles come from turgor pressure P , and the tension j p . The shape of the cell poles are then 
described by the two-dimensional Laplace’s law, 


7 f' sw c±{u,t) = P , (S.25) 

where the superscripts st and suu denote respectively the stalked and swarmer poles. 

(2) Upper cell contour (l s t/ sm ): The region l st / sw in the upper cell contour obeys the force-balance equation: 


P 7 C + (+ f) T /width(n, t) + /cres(u, t) — 0 . 


(S.26) 
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At the endpoints of the segments l s t/ S w we impose the boundary condition that the curvature c must equal the 
longitudinal curvature of the poles. 

(3) Bottom cell contour (2 st / sw ): In the absence of crescentin, the region 2 st / sw in the bottom cell contour obeys 
the force-balance equation: 


P ~ 7 C - (u,t) + /width (u,t) = 0 . (S.27) 

(4) Septal region: We incorporate the effect of constriction in the cell shape equation by imposing the constraint 
(boundary condition) that w(l st (t),t) = w m i n (t), where w m i n (t) evolves according to the kinetics described in Eq. ( 6 ) 
and (7) of the main text. Furthermore, the curvatures at the end points of the septal segments must conform to the 
curvature of the newly formed poles. Subject to these boundary conditions, the width profile in the upper septal 
region ( I st - l a /2 < u < I st + l a /2) is given by, 


P - 7 s C + {u,t) + /width(w> t) + fcres{u,t) = 0 , (S.28) 

with a surface tension (y s ) of the newly formed poles chosen to be much higher than the peripheral region, ~ 4.5q. 
The contour below the centerline obeys the equation, 


p - 7 sC-(u, t) + /width (u, t) = 0 . 


(S.29) 


The governing cell shape equation given in Eq. ( |S.20 ), is then solved numerically in each part of the cell contour 
with matching boundary conditions in w±{u,t ) and its derivatives. 
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